This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.2419355

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.12 while 0.88 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

11

1

1

0.474

0.468

0.237

0.728

1.026

0.127

0.047

CSEO

p_pelagic

11

1

2

0.409

0.403

0.192

0.651

1.011

0.122

0.049

CSEO

p_pelagic

11

10

1

0.477

0.476

0.233

0.724

1.022

0.126

0.008

CSEO

p_pelagic

11

10

2

0.412

0.403

0.189

0.668

1.022

0.124

0.073

CSEO

p_pelagic

11

11

1

0.471

0.468

0.242

0.713

1.021

0.122

0.025

CSEO

p_pelagic

11

11

2

0.409

0.402

0.191

0.668

1.022

0.125

0.056

CSEO

p_pelagic

11

12

1

0.475

0.472

0.237

0.717

1.015

0.127

0.024

CSEO

p_pelagic

11

12

2

0.406

0.398

0.192

0.650

1.017

0.121

0.066

CSEO

p_yellow

11

12

1

0.518

0.517

0.235

0.802

1.011

0.149

0.007

CSEO

p_pelagic

11

13

1

0.476

0.476

0.232

0.730

1.014

0.128

0.000

CSEO

p_pelagic

11

13

2

0.409

0.404

0.187

0.655

1.015

0.123

0.041

CSEO

p_pelagic

11

14

1

0.485

0.485

0.252

0.730

1.028

0.125

0.000

CSEO

p_pelagic

11

14

2

0.414

0.408

0.193

0.671

1.018

0.123

0.049

CSEO

p_pelagic

11

15

1

0.444

0.440

0.216

0.705

1.030

0.127

0.032

CSEO

p_pelagic

11

15

2

0.403

0.396

0.187

0.669

1.015

0.124

0.056

CSEO

p_pelagic

11

16

1

0.491

0.489

0.258

0.723

1.022

0.122

0.016

CSEO

p_pelagic

11

16

2

0.414

0.409

0.185

0.665

1.025

0.125

0.040

CSEO

p_pelagic

11

17

1

0.481

0.477

0.256

0.714

1.011

0.121

0.033

CSEO

p_pelagic

11

17

2

0.410

0.404

0.187

0.669

1.022

0.123

0.049

CSEO

p_pelagic

11

18

1

0.481

0.480

0.245

0.726

1.019

0.123

0.008

CSEO

p_pelagic

11

18

2

0.411

0.410

0.187

0.658

1.010

0.122

0.008

CSEO

p_pelagic

11

19

1

0.494

0.494

0.262

0.728

1.022

0.121

0.000

CSEO

p_pelagic

11

19

2

0.416

0.408

0.190

0.675

1.018

0.126

0.063

CSEO

p_pelagic

11

2

1

0.474

0.473

0.242

0.719

1.019

0.127

0.008

CSEO

p_pelagic

11

2

2

0.410

0.404

0.188

0.664

1.011

0.125

0.048

CSEO

p_yellow

11

2

2

0.426

0.413

0.169

0.737

1.011

0.150

0.086

CSEO

p_pelagic

11

20

1

0.488

0.483

0.263

0.733

1.029

0.123

0.041

CSEO

p_pelagic

11

20

2

0.416

0.411

0.193

0.672

1.017

0.124

0.040

CSEO

p_pelagic

11

21

2

0.415

0.409

0.189

0.678

1.011

0.125

0.048

CSEO

p_yellow

11

21

2

0.428

0.419

0.162

0.750

1.013

0.154

0.058

CSEO

p_pelagic

11

22

2

0.417

0.411

0.196

0.673

1.016

0.123

0.049

CSEO

p_yellow

11

22

1

0.521

0.519

0.242

0.809

1.014

0.151

0.013

CSEO

p_pelagic

11

23

2

0.452

0.449

0.233

0.696

1.012

0.121

0.025

CSEO

p_yellow

11

23

1

0.537

0.536

0.268

0.812

1.011

0.143

0.007

CSEO

p_yellow

11

25

1

0.486

0.478

0.222

0.776

1.011

0.145

0.055

CSEO

p_yellow

11

25

2

0.442

0.432

0.170

0.758

1.011

0.154

0.065

CSEO

p_yellow

11

26

1

0.512

0.509

0.243

0.794

1.011

0.145

0.021

CSEO

p_yellow

11

28

1

0.486

0.481

0.214

0.774

1.011

0.149

0.033

CSEO

p_yellow

11

29

1

0.501

0.499

0.230

0.800

1.014

0.148

0.014

CSEO

p_pelagic

11

3

1

0.475

0.471

0.237

0.724

1.021

0.128

0.031

CSEO

p_pelagic

11

3

2

0.409

0.403

0.182

0.671

1.013

0.126

0.048

CSEO

p_pelagic

11

4

1

0.479

0.478

0.248

0.723

1.013

0.125

0.008

CSEO

p_pelagic

11

4

2

0.407

0.401

0.187

0.650

1.017

0.121

0.050

CSEO

p_yellow

11

4

2

0.428

0.418

0.175

0.738

1.016

0.147

0.068

CSEO

p_pelagic

11

5

1

0.476

0.474

0.238

0.720

1.011

0.127

0.016

CSEO

p_pelagic

11

5

2

0.411

0.402

0.189

0.674

1.021

0.124

0.073

CSEO

p_pelagic

11

6

1

0.481

0.479

0.247

0.730

1.014

0.126

0.016

CSEO

p_pelagic

11

6

2

0.412

0.408

0.185

0.667

1.023

0.124

0.032

CSEO

p_yellow

11

6

1

0.516

0.516

0.238

0.806

1.011

0.151

0.000

CSEO

p_pelagic

11

7

1

0.476

0.475

0.233

0.722

1.022

0.126

0.008

CSEO

p_pelagic

11

7

2

0.411

0.407

0.192

0.666

1.023

0.123

0.033

CSEO

p_pelagic

11

8

1

0.478

0.479

0.241

0.731

1.015

0.128

-0.008

CSEO

p_pelagic

11

8

2

0.407

0.401

0.189

0.654

1.019

0.123

0.049

CSEO

p_pelagic

11

9

1

0.474

0.470

0.237

0.727

1.017

0.128

0.031

CSEO

p_pelagic

11

9

2

0.411

0.409

0.189

0.662

1.024

0.125

0.016

CSEO

p_pelagic

13

1

1

0.343

0.328

0.123

0.634

1.074

0.135

0.111

NSEI

p_pelagic

13

1

2

0.149

0.132

0.043

0.346

1.066

0.079

0.216

NSEI

p_yellow

13

1

1

0.430

0.374

0.129

0.982

1.029

0.224

0.250

NSEI

p_yellow

13

1

2

0.431

0.380

0.127

0.981

1.028

0.223

0.228

NSEI

p_pelagic

13

10

1

0.352

0.341

0.125

0.637

1.060

0.132

0.083

NSEI

p_pelagic

13

10

2

0.155

0.138

0.046

0.357

1.074

0.081

0.211

NSEI

p_yellow

13

10

1

0.429

0.381

0.125

0.982

1.032

0.224

0.214

NSEI

p_yellow

13

10

2

0.433

0.384

0.126

0.984

1.029

0.224

0.219

NSEI

p_pelagic

13

11

1

0.354

0.340

0.134

0.647

1.053

0.133

0.105

NSEI

p_pelagic

13

11

2

0.154

0.138

0.046

0.348

1.049

0.079

0.202

NSEI

p_yellow

13

11

1

0.433

0.380

0.126

0.982

1.028

0.225

0.236

NSEI

p_yellow

13

11

2

0.434

0.384

0.129

0.983

1.031

0.224

0.224

NSEI

p_pelagic

13

12

1

0.360

0.350

0.140

0.632

1.048

0.131

0.076

NSEI

p_pelagic

13

12

2

0.154

0.141

0.045

0.340

1.064

0.078

0.167

NSEI

p_yellow

13

12

1

0.431

0.381

0.127

0.983

1.030

0.224

0.223

NSEI

p_yellow

13

12

2

0.434

0.380

0.132

0.981

1.030

0.226

0.239

NSEI

p_pelagic

13

13

1

0.365

0.357

0.136

0.642

1.057

0.133

0.060

NSEI

p_pelagic

13

13

2

0.158

0.146

0.048

0.346

1.046

0.079

0.153

NSEI

p_yellow

13

13

1

0.433

0.384

0.135

0.980

1.029

0.223

0.219

NSEI

p_yellow

13

13

2

0.432

0.379

0.130

0.985

1.031

0.224

0.236

NSEI

p_pelagic

13

14

1

0.343

0.328

0.139

0.621

1.052

0.124

0.121

NSEI

p_pelagic

13

14

2

0.155

0.141

0.049

0.340

1.076

0.076

0.185

NSEI

p_yellow

13

14

1

0.428

0.370

0.131

0.981

1.028

0.225

0.258

NSEI

p_yellow

13

14

2

0.436

0.384

0.129

0.981

1.027

0.225

0.232

NSEI

p_pelagic

13

15

1

0.315

0.299

0.120

0.578

1.044

0.120

0.133

NSEI

p_pelagic

13

15

2

0.145

0.133

0.046

0.321

1.063

0.072

0.166

NSEI

p_yellow

13

15

1

0.426

0.376

0.124

0.982

1.034

0.224

0.223

NSEI

p_yellow

13

15

2

0.432

0.380

0.131

0.985

1.032

0.224

0.233

NSEI

p_pelagic

13

16

1

0.350

0.341

0.141

0.616

1.038

0.125

0.072

NSEI

p_pelagic

13

16

2

0.162

0.149

0.053

0.350

1.046

0.078

0.166

NSEI

p_yellow

13

16

1

0.429

0.377

0.128

0.979

1.032

0.224

0.233

NSEI

p_yellow

13

16

2

0.431

0.383

0.127

0.983

1.026

0.224

0.215

NSEI

p_pelagic

13

17

1

0.363

0.355

0.159

0.620

1.030

0.123

0.065

NSEI

p_pelagic

13

17

2

0.163

0.151

0.055

0.340

1.062

0.076

0.159

NSEI

p_yellow

13

17

1

0.431

0.381

0.126

0.982

1.033

0.225

0.222

NSEI

p_yellow

13

17

2

0.434

0.382

0.132

0.985

1.029

0.224

0.232

NSEI

p_pelagic

13

18

1

0.368

0.356

0.159

0.631

1.041

0.124

0.097

NSEI

p_pelagic

13

18

2

0.167

0.154

0.056

0.349

1.038

0.076

0.170

NSEI

p_yellow

13

18

1

0.432

0.382

0.121

0.981

1.030

0.224

0.223

NSEI

p_yellow

13

18

2

0.434

0.381

0.132

0.981

1.030

0.224

0.237

NSEI

p_pelagic

13

19

1

0.366

0.353

0.156

0.633

1.036

0.123

0.106

NSEI

p_pelagic

13

19

2

0.170

0.158

0.059

0.355

1.035

0.078

0.153

NSEI

p_yellow

13

19

1

0.428

0.381

0.124

0.980

1.031

0.223

0.211

NSEI

p_yellow

13

19

2

0.434

0.385

0.134

0.982

1.028

0.225

0.218

NSEI

p_pelagic

13

2

1

0.343

0.328

0.129

0.628

1.077

0.133

0.113

NSEI

p_pelagic

13

2

2

0.149

0.132

0.043

0.337

1.083

0.078

0.217

NSEI

p_yellow

13

2

1

0.431

0.382

0.127

0.982

1.033

0.224

0.219

NSEI

p_yellow

13

2

2

0.435

0.386

0.130

0.981

1.028

0.224

0.219

NSEI

p_pelagic

13

20

1

0.380

0.372

0.167

0.651

1.032

0.123

0.065

NSEI

p_pelagic

13

20

2

0.176

0.162

0.061

0.353

1.027

0.076

0.185

NSEI

p_yellow

13

20

1

0.426

0.376

0.124

0.981

1.031

0.224

0.223

NSEI

p_yellow

13

20

2

0.430

0.375

0.126

0.981

1.025

0.225

0.244

NSEI

p_pelagic

13

21

1

0.397

0.390

0.186

0.649

1.031

0.122

0.057

NSEI

p_pelagic

13

21

2

0.185

0.171

0.065

0.377

1.023

0.082

0.170

NSEI

p_yellow

13

21

1

0.428

0.378

0.133

0.983

1.035

0.223

0.225

NSEI

p_yellow

13

21

2

0.431

0.381

0.131

0.984

1.027

0.224

0.223

NSEI

p_pelagic

13

22

2

0.184

0.173

0.064

0.365

1.016

0.078

0.141

NSEI

p_yellow

13

22

1

0.424

0.373

0.122

0.980

1.038

0.224

0.228

NSEI

p_yellow

13

22

2

0.432

0.381

0.128

0.983

1.025

0.223

0.228

NSEI

p_pelagic

13

23

2

0.117

0.108

0.043

0.242

1.016

0.051

0.175

NSEI

p_yellow

13

23

1

0.392

0.335

0.103

0.974

1.030

0.225

0.254

NSEI

p_yellow

13

23

2

0.483

0.446

0.165

0.987

1.018

0.216

0.171

NSEI

p_yellow

13

24

1

0.395

0.339

0.103

0.976

1.035

0.228

0.246

NSEI

p_yellow

13

24

2

0.481

0.441

0.162

0.986

1.021

0.217

0.184

NSEI

p_yellow

13

25

1

0.396

0.341

0.108

0.975

1.031

0.223

0.246

NSEI

p_yellow

13

25

2

0.480

0.441

0.154

0.986

1.024

0.217

0.180

NSEI

p_yellow

13

26

1

0.387

0.330

0.098

0.975

1.027

0.224

0.255

NSEI

p_yellow

13

26

2

0.471

0.433

0.158

0.986

1.027

0.218

0.175

NSEI

p_yellow

13

27

1

0.380

0.326

0.098

0.967

1.026

0.219

0.247

NSEI

p_yellow

13

27

2

0.473

0.436

0.155

0.984

1.023

0.214

0.173

NSEI

p_yellow

13

28

1

0.388

0.342

0.101

0.968

1.025

0.215

0.214

NSEI

p_yellow

13

28

2

0.476

0.440

0.161

0.981

1.029

0.208

0.173

NSEI

p_yellow

13

29

1

0.365

0.318

0.099

0.936

1.024

0.200

0.235

NSEI

p_yellow

13

29

2

0.468

0.442

0.156

0.956

1.020

0.200

0.130

NSEI

p_pelagic

13

3

1

0.344

0.329

0.130

0.634

1.074

0.133

0.113

NSEI

p_pelagic

13

3

2

0.150

0.133

0.042

0.343

1.075

0.079

0.214

NSEI

p_yellow

13

3

1

0.431

0.378

0.127

0.980

1.031

0.224

0.237

NSEI

p_yellow

13

3

2

0.434

0.385

0.130

0.981

1.029

0.225

0.218

NSEI

p_pelagic

13

4

1

0.344

0.334

0.129

0.622

1.054

0.131

0.076

NSEI

p_pelagic

13

4

2

0.149

0.133

0.042

0.349

1.084

0.079

0.203

NSEI

p_yellow

13

4

1

0.433

0.386

0.128

0.983

1.027

0.224

0.209

NSEI

p_yellow

13

4

2

0.439

0.385

0.134

0.984

1.033

0.225

0.240

NSEI

p_pelagic

13

5

1

0.345

0.332

0.125

0.622

1.053

0.131

0.099

NSEI

p_pelagic

13

5

2

0.150

0.134

0.042

0.342

1.054

0.079

0.202

NSEI

p_yellow

13

5

1

0.429

0.378

0.129

0.984

1.032

0.224

0.228

NSEI

p_yellow

13

5

2

0.434

0.380

0.136

0.984

1.033

0.225

0.240

NSEI

p_pelagic

13

6

1

0.350

0.337

0.131

0.627

1.057

0.133

0.098

NSEI

p_pelagic

13

6

2

0.150

0.134

0.042

0.352

1.077

0.079

0.202

NSEI

p_yellow

13

6

1

0.430

0.378

0.128

0.980

1.035

0.226

0.230

NSEI

p_yellow

13

6

2

0.432

0.383

0.131

0.984

1.030

0.225

0.218

NSEI

p_pelagic

13

7

1

0.350

0.340

0.130

0.631

1.053

0.133

0.075

NSEI

p_pelagic

13

7

2

0.153

0.134

0.045

0.345

1.058

0.081

0.236

NSEI

p_yellow

13

7

1

0.431

0.385

0.128

0.982

1.027

0.224

0.206

NSEI

p_yellow

13

7

2

0.438

0.389

0.131

0.982

1.030

0.224

0.218

NSEI

p_pelagic

13

8

1

0.350

0.339

0.131

0.625

1.062

0.132

0.083

NSEI

p_pelagic

13

8

2

0.151

0.135

0.044

0.348

1.072

0.079

0.203

NSEI

p_yellow

13

8

1

0.429

0.379

0.121

0.982

1.031

0.223

0.224

NSEI

p_yellow

13

8

2

0.441

0.395

0.135

0.984

1.033

0.224

0.205

NSEI

p_pelagic

13

9

1

0.349

0.335

0.133

0.631

1.075

0.129

0.108

NSEI

p_pelagic

13

9

2

0.152

0.138

0.045

0.342

1.081

0.077

0.181

NSEI

p_yellow

13

9

1

0.434

0.382

0.130

0.984

1.030

0.224

0.232

NSEI

p_yellow

13

9

2

0.435

0.389

0.127

0.982

1.029

0.224

0.205

NSEI

p_yellow

14

1

1

0.585

0.567

0.256

0.985

1.010

0.189

0.095

NSEO

p_yellow

14

10

1

0.579

0.563

0.254

0.985

1.015

0.189

0.085

NSEO

p_yellow

14

11

1

0.587

0.572

0.254

0.985

1.010

0.191

0.079

NSEO

p_yellow

14

12

1

0.587

0.572

0.266

0.984

1.013

0.189

0.079

NSEO

p_yellow

14

12

2

0.622

0.614

0.287

0.989

1.011

0.185

0.043

NSEO

p_yellow

14

13

1

0.586

0.569

0.251

0.984

1.014

0.188

0.090

NSEO

p_yellow

14

14

1

0.586

0.572

0.258

0.988

1.012

0.190

0.074

NSEO

p_yellow

14

14

2

0.631

0.622

0.303

0.989

1.011

0.180

0.050

NSEO

p_yellow

14

16

2

0.627

0.621

0.299

0.988

1.011

0.182

0.033

NSEO

p_yellow

14

18

1

0.583

0.564

0.252

0.986

1.011

0.192

0.099

NSEO

p_yellow

14

18

2

0.624

0.614

0.285

0.989

1.011

0.184

0.054

NSEO

p_yellow

14

19

1

0.589

0.575

0.252

0.986

1.010

0.191

0.073

NSEO

p_yellow

14

19

2

0.625

0.617

0.288

0.989

1.014

0.181

0.044

NSEO

p_yellow

14

2

1

0.584

0.567

0.262

0.985

1.013

0.188

0.090

NSEO

p_yellow

14

20

1

0.586

0.572

0.251

0.985

1.013

0.190

0.074

NSEO

p_yellow

14

21

1

0.589

0.578

0.255

0.987

1.014

0.189

0.058

NSEO

p_yellow

14

22

2

0.625

0.618

0.296

0.991

1.012

0.181

0.039

NSEO

p_yellow

14

23

1

0.618

0.611

0.299

0.986

1.013

0.178

0.039

NSEO

p_yellow

14

23

2

0.628

0.620

0.284

0.988

1.010

0.183

0.044

NSEO

p_yellow

14

26

1

0.595

0.585

0.264

0.985

1.012

0.188

0.053

NSEO

p_yellow

14

26

2

0.625

0.620

0.287

0.990

1.012

0.184

0.027

NSEO

p_yellow

14

28

1

0.591

0.581

0.272

0.985

1.017

0.182

0.055

NSEO

p_yellow

14

28

2

0.621

0.614

0.284

0.988

1.010

0.184

0.038

NSEO

p_yellow

14

29

1

0.592

0.588

0.262

0.976

1.012

0.180

0.022

NSEO

p_yellow

14

29

2

0.618

0.615

0.282

0.982

1.013

0.180

0.017

NSEO

p_yellow

14

4

1

0.586

0.571

0.254

0.984

1.012

0.190

0.079

NSEO

p_dsr

14

44

1

0.076

0.059

0.015

0.234

1.014

0.058

0.293

NSEO

p_yellow

14

5

2

0.623

0.613

0.274

0.989

1.010

0.185

0.054

NSEO

p_yellow

14

6

1

0.588

0.572

0.264

0.986

1.012

0.189

0.085

NSEO

p_yellow

14

6

2

0.629

0.621

0.306

0.987

1.013

0.180

0.044

NSEO

p_yellow

14

8

1

0.586

0.570

0.254

0.985

1.013

0.189

0.085

NSEO

p_pelagic

16

1

1

0.380

0.368

0.140

0.664

1.062

0.137

0.088

SSEO

p_pelagic

16

1

2

0.359

0.348

0.125

0.629

1.052

0.134

0.082

SSEO

p_pelagic

16

10

1

0.386

0.377

0.157

0.671

1.045

0.135

0.067

SSEO

p_pelagic

16

10

2

0.366

0.357

0.140

0.652

1.040

0.134

0.067

SSEO

p_pelagic

16

11

1

0.386

0.378

0.156

0.663

1.049

0.134

0.060

SSEO

p_pelagic

16

11

2

0.366

0.362

0.140

0.637

1.038

0.133

0.030

SSEO

p_pelagic

16

12

1

0.390

0.379

0.153

0.668

1.044

0.136

0.081

SSEO

p_pelagic

16

12

2

0.364

0.358

0.139

0.636

1.041

0.131

0.046

SSEO

p_pelagic

16

13

1

0.392

0.384

0.156

0.657

1.042

0.134

0.060

SSEO

p_pelagic

16

13

2

0.367

0.358

0.144

0.641

1.039

0.133

0.068

SSEO

p_pelagic

16

14

1

0.375

0.366

0.156

0.651

1.032

0.130

0.069

SSEO

p_pelagic

16

14

2

0.366

0.357

0.142

0.653

1.043

0.132

0.068

SSEO

p_pelagic

16

15

1

0.381

0.374

0.162

0.646

1.029

0.128

0.055

SSEO

p_pelagic

16

15

2

0.369

0.360

0.144

0.635

1.030

0.129

0.070

SSEO

p_pelagic

16

16

1

0.382

0.373

0.153

0.645

1.033

0.131

0.069

SSEO

p_pelagic

16

16

2

0.370

0.360

0.144

0.636

1.042

0.128

0.078

SSEO

p_pelagic

16

17

1

0.427

0.421

0.194

0.688

1.042

0.127

0.047

SSEO

p_pelagic

16

17

2

0.383

0.375

0.154

0.658

1.024

0.129

0.062

SSEO

p_pelagic

16

18

1

0.416

0.410

0.187

0.668

1.028

0.127

0.047

SSEO

p_pelagic

16

18

2

0.386

0.377

0.161

0.648

1.024

0.128

0.070

SSEO

p_pelagic

16

19

1

0.399

0.392

0.185

0.651

1.021

0.123

0.057

SSEO

p_pelagic

16

19

2

0.385

0.379

0.163

0.647

1.027

0.128

0.047

SSEO

p_pelagic

16

2

1

0.381

0.374

0.137

0.664

1.062

0.139

0.050

SSEO

p_pelagic

16

2

2

0.362

0.352

0.132

0.646

1.056

0.136

0.074

SSEO

p_pelagic

16

20

1

0.404

0.396

0.185

0.658

1.016

0.123

0.065

SSEO

p_pelagic

16

20

2

0.388

0.379

0.166

0.648

1.012

0.128

0.070

SSEO

p_pelagic

16

21

1

0.388

0.379

0.172

0.642

1.013

0.123

0.073

SSEO

p_pelagic

16

21

2

0.385

0.379

0.169

0.648

1.013

0.125

0.048

SSEO

p_pelagic

16

3

1

0.379

0.370

0.140

0.662

1.052

0.137

0.066

SSEO

p_pelagic

16

3

2

0.362

0.352

0.130

0.653

1.056

0.136

0.073

SSEO

p_dsr

16

31

1

0.402

0.400

0.375

0.434

1.783

0.016

0.129

SSEO

p_dsr

16

36

1

0.603

0.603

0.573

0.634

1.014

0.015

0.000

SSEO

p_slope

16

36

1

0.060

0.061

0.040

0.078

1.068

0.009

-0.107

SSEO

p_dsr

16

37

1

0.616

0.616

0.587

0.646

1.013

0.015

0.000

SSEO

p_slope

16

37

1

0.065

0.065

0.044

0.083

1.089

0.010

0.000

SSEO

p_pelagic

16

4

1

0.381

0.370

0.144

0.659

1.047

0.135

0.081

SSEO

p_pelagic

16

4

2

0.359

0.354

0.130

0.624

1.051

0.132

0.038

SSEO

p_pelagic

16

42

1

0.833

0.833

0.826

0.840

1.111

0.003

0.000

SSEO

p_yellow

16

42

1

0.407

0.407

0.385

0.428

1.038

0.011

0.000

SSEO

p_dsr

16

42

1

0.451

0.441

0.405

0.509

3.601

0.031

0.326

SSEO

p_slope

16

42

1

0.167

0.180

0.109

0.210

4.888

0.032

-0.405

SSEO

p_pelagic

16

5

1

0.382

0.375

0.148

0.668

1.050

0.138

0.051

SSEO

p_pelagic

16

5

2

0.357

0.349

0.127

0.633

1.059

0.133

0.060

SSEO

p_pelagic

16

6

1

0.385

0.378

0.141

0.658

1.058

0.136

0.052

SSEO

p_pelagic

16

6

2

0.361

0.355

0.133

0.635

1.059

0.133

0.045

SSEO

p_pelagic

16

7

1

0.378

0.375

0.138

0.657

1.047

0.135

0.022

SSEO

p_pelagic

16

7

2

0.362

0.349

0.132

0.637

1.066

0.135

0.096

SSEO

p_pelagic

16

8

1

0.388

0.379

0.146

0.662

1.052

0.135

0.066

SSEO

p_pelagic

16

8

2

0.361

0.354

0.133

0.646

1.057

0.134

0.052

SSEO

p_pelagic

16

9

1

0.384

0.378

0.151

0.657

1.046

0.134

0.045

SSEO

p_pelagic

16

9

2

0.362

0.353

0.136

0.631

1.042

0.133

0.068

SSEO

p_pelagic

4

1

1

0.559

0.564

0.292

0.793

1.029

0.130

-0.039

PWSO

p_pelagic

4

1

2

0.555

0.565

0.270

0.798

1.032

0.134

-0.075

PWSO

p_pelagic

4

10

1

0.565

0.573

0.300

0.794

1.029

0.130

-0.062

PWSO

p_pelagic

4

10

2

0.559

0.565

0.289

0.792

1.030

0.130

-0.046

PWSO

p_pelagic

4

11

1

0.567

0.574

0.309

0.799

1.038

0.129

-0.054

PWSO

p_pelagic

4

11

2

0.557

0.562

0.294

0.793

1.028

0.131

-0.038

PWSO

p_pelagic

4

12

1

0.558

0.558

0.295

0.793

1.030

0.132

0.000

PWSO

p_pelagic

4

12

2

0.554

0.557

0.285

0.795

1.025

0.133

-0.023

PWSO

p_pelagic

4

13

1

0.566

0.574

0.295

0.809

1.034

0.132

-0.060

PWSO

p_pelagic

4

13

2

0.560

0.567

0.295

0.791

1.028

0.130

-0.054

PWSO

p_pelagic

4

14

1

0.561

0.566

0.302

0.795

1.029

0.130

-0.038

PWSO

p_pelagic

4

14

2

0.557

0.562

0.297

0.786

1.023

0.131

-0.038

PWSO

p_pelagic

4

15

1

0.559

0.562

0.292

0.797

1.033

0.131

-0.023

PWSO

p_pelagic

4

15

2

0.561

0.565

0.305

0.804

1.027

0.128

-0.031

PWSO

p_pelagic

4

16

1

0.559

0.564

0.292

0.783

1.021

0.128

-0.039

PWSO

p_pelagic

4

16

2

0.555

0.560

0.289

0.795

1.031

0.130

-0.038

PWSO

p_pelagic

4

17

1

0.562

0.565

0.299

0.800

1.038

0.131

-0.023

PWSO

p_pelagic

4

17

2

0.558

0.564

0.299

0.796

1.027

0.131

-0.046

PWSO

p_pelagic

4

18

1

0.561

0.564

0.300

0.798

1.028

0.130

-0.023

PWSO

p_pelagic

4

18

2

0.558

0.564

0.291

0.800

1.019

0.132

-0.045

PWSO

p_pelagic

4

19

1

0.555

0.558

0.290

0.798

1.027

0.131

-0.023

PWSO

p_pelagic

4

19

2

0.557

0.567

0.300

0.797

1.038

0.132

-0.076

PWSO

p_pelagic

4

2

1

0.562

0.570

0.298

0.792

1.032

0.131

-0.061

PWSO

p_pelagic

4

2

2

0.557

0.562

0.291

0.794

1.030

0.131

-0.038

PWSO

p_pelagic

4

20

2

0.471

0.468

0.244

0.708

1.023

0.121

0.025

PWSO

p_pelagic

4

21

2

0.410

0.409

0.245

0.589

1.017

0.088

0.011

PWSO

p_pelagic

4

22

2

0.465

0.465

0.284

0.652

1.015

0.095

0.000

PWSO

p_pelagic

4

3

1

0.565

0.575

0.285

0.803

1.035

0.132

-0.076

PWSO

p_pelagic

4

3

2

0.555

0.563

0.292

0.794

1.035

0.132

-0.061

PWSO

p_pelagic

4

4

1

0.562

0.568

0.296

0.798

1.017

0.129

-0.046

PWSO

p_pelagic

4

4

2

0.555

0.562

0.287

0.803

1.038

0.134

-0.052

PWSO

p_pelagic

4

5

1

0.562

0.570

0.307

0.792

1.032

0.129

-0.062

PWSO

p_pelagic

4

5

2

0.555

0.559

0.294

0.787

1.022

0.128

-0.031

PWSO

p_pelagic

4

6

1

0.563

0.569

0.307

0.790

1.032

0.128

-0.047

PWSO

p_pelagic

4

6

2

0.560

0.566

0.299

0.795

1.029

0.130

-0.046

PWSO

p_pelagic

4

7

1

0.562

0.571

0.295

0.795

1.022

0.132

-0.068

PWSO

p_pelagic

4

7

2

0.560

0.565

0.294

0.793

1.033

0.130

-0.039

PWSO

p_pelagic

4

8

1

0.561

0.567

0.301

0.795

1.031

0.129

-0.046

PWSO

p_pelagic

4

8

2

0.556

0.561

0.291

0.798

1.038

0.133

-0.037

PWSO

p_pelagic

4

9

1

0.567

0.572

0.297

0.800

1.027

0.129

-0.039

PWSO

p_pelagic

4

9

2

0.560

0.565

0.298

0.801

1.041

0.133

-0.038

PWSO

p_black

8

40

1

0.440

0.412

0.340

0.949

1.126

0.128

0.219

afognak

p_black

8

48

1

0.525

0.491

0.414

0.963

1.053

0.129

0.263

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta0_pH

10

1

0.265

0.292

-0.231

0.616

1.119

0.208

-0.130

northeast

beta1_pH

10

1

2.174

2.058

1.460

4.570

1.258

0.615

0.189

northeast

beta2_pH

10

1

1.706

0.912

0.140

9.578

1.319

2.330

0.341

northeast

beta2_pH

10

3

0.868

0.277

0.105

11.247

1.243

2.335

0.253

northeast

beta3_pH

10

1

32.532

32.273

30.069

38.076

1.182

1.738

0.149

northeast

beta0_pH

11

1

0.103

0.108

-0.551

0.676

1.132

0.310

-0.016

CSEO

beta2_pH

11

1

0.467

0.408

0.221

1.053

1.104

0.226

0.260

CSEO

beta0_pelagic

11

-0.101

-0.093

-0.473

0.196

1.284

0.166

-0.048

CSEO

beta1_pelagic

11

1.088

1.090

0.678

1.531

1.282

0.227

-0.009

CSEO

beta2_pelagic

11

1.870

1.377

0.303

6.032

1.179

1.534

0.321

CSEO

beta0_pelagic

13

-0.728

-0.698

-1.442

-0.019

1.226

0.389

-0.077

NSEI

beta1_pelagic

13

2.368

2.340

1.346

3.395

1.215

0.563

0.050

NSEI

beta3_pelagic

13

31.557

31.037

25.510

39.680

1.184

3.725

0.140

NSEI

beta3_pelagic

14

34.340

34.189

31.378

38.454

1.102

1.739

0.087

NSEO

beta0_pH

15

1

0.086

0.103

-0.333

0.405

1.108

0.186

-0.091

SSEI

beta0_yellow

15

-1.261

-1.235

-1.785

-1.022

1.175

0.174

-0.149

SSEI

beta1_yellow

15

0.175

0.079

0.011

0.906

1.193

0.232

0.413

SSEI

beta2_yellow

15

-5.846

-4.131

-20.009

-0.182

1.330

5.431

-0.316

SSEI

beta0_pH

16

5

0.252

0.240

0.157

0.412

1.136

0.067

0.178

SSEO

beta3_pH

16

4

41.579

41.585

40.888

42.388

1.136

0.402

-0.015

SSEO

beta0_pelagic

16

-0.539

-0.485

-1.399

-0.004

1.230

0.383

-0.141

SSEO

beta1_pelagic

16

1.646

1.510

0.816

3.002

1.175

0.599

0.227

SSEO

beta0_dsr

16

0.160

0.160

0.128

0.193

1.224

0.016

0.000

SSEO

beta0_slope

16

-2.650

-2.651

-2.730

-2.562

1.174

0.043

0.024

SSEO

mu_beta2_pH

2

1

1.762

0.941

0.163

9.219

1.377

2.490

0.330

NG

beta2_yellow

2

-9.401

-9.118

-20.733

-0.522

1.216

5.790

-0.049

NG

beta1_pH

4

1

1.960

1.882

1.246

3.700

1.261

0.535

0.146

PWSO

beta2_pH

4

1

1.029

0.790

0.194

2.896

1.185

0.741

0.323

PWSO

beta0_pelagic

4

0.268

0.293

-0.292

0.702

1.214

0.243

-0.103

PWSO

beta1_pelagic

4

0.778

0.755

0.267

1.379

1.202

0.274

0.084

PWSO

beta0_yellow

4

0.844

0.949

-0.124

1.188

1.138

0.333

-0.315

PWSO

beta2_pH

5

1

2.357

1.570

0.282

10.127

1.346

2.662

0.296

BSAI

beta0_pH

6

1

-0.449

-0.329

-1.603

0.030

1.253

0.405

-0.296

SOKO2SAP

beta1_pH

6

1

1.933

1.788

1.080

3.367

1.101

0.621

0.234

SOKO2SAP

beta2_pH

6

1

1.644

0.577

0.141

10.824

1.367

2.816

0.379

SOKO2SAP

beta3_pH

6

1

33.025

33.148

25.836

38.926

1.137

2.871

-0.043

SOKO2SAP

beta1_pH

8

1

2.085

1.944

1.076

3.838

1.151

0.703

0.200

afognak

beta2_pH

8

1

1.363

0.465

0.150

9.360

1.436

2.309

0.389

afognak

beta2_pH

8

3

1.976

1.135

0.143

11.512

1.174

2.524

0.333

afognak

beta1_pH

9

1

1.963

1.887

1.240

3.078

1.147

0.471

0.161

eastside

beta2_pH

9

1

1.437

0.621

0.188

8.572

1.352

2.235

0.365

eastside

beta2_pH

9

3

0.761

0.213

0.092

11.002

1.269

2.335

0.235

eastside

lambda_H

5

2.908

0.621

0.026

21.835

1.139

8.713

0.262

BSAI

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

811.00

274.00

1,215.00

347.00

Prop. of total

0.75

0.25

0.78

0.22

Prop. median less than 1

0.17

0.19

0.28

0.41

Prop. median less than 10

0.24

0.47

0.44

0.70

Prop. median less than 50

0.34

0.69

0.60

0.78

Median of the median values

402.27

12.24

18.77

1.90

Mean of the median values

1,437.40

561.68

452.16

187.14

Mean skew of posteriors

0.20

0.28

0.26

0.26

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

12.00

11.00

59.00

38.00

Prop. of total

0.52

0.48

0.61

0.39

Prop. median less than 1

0.25

0.55

0.64

0.66

Prop. median less than 10

0.25

1.00

0.69

0.84

Prop. median less than 50

0.25

1.00

0.83

0.89

Median of the median values

286.28

0.98

0.01

0.09

Mean of the median values

635.64

1.12

87.35

19.43

Mean skew of posteriors

0.22

0.21

0.15

0.18

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_dsr

1

beta0_pelagic

4

beta0_pH

4

beta0_slope

1

beta0_yellow

2

beta1_pelagic

4

beta1_pH

4

beta1_yellow

1

beta2_pelagic

1

beta2_pH

9

beta2_yellow

2

beta3_pelagic

1

beta3_pH

3

H_ay

3

H_ayg

1

H_ayu

3

Hb_ay

3

Hb_ayg

2

Hb_ayu

3

Hd_ay

1

Hd_ayg

2

Ho_ay

1

Ho_ayg

1

Ho_ayu

3

Hp_ay

3

Hp_ayg

1

Hp_ayu

3

Hs_ay

2

Hs_ayg

2

Hs_ayu

1

Hy_ay

4

Hy_ayg

2

Hy_ayu

5

lambda_H

1

mu_beta2_pH

1

p_black

1

p_dsr

2

p_pelagic

1

p_slope

2

pH

16

R_ay

3

R_ayg

3

R_ayu

5

Rb_ay

3

Rb_ayg

4

Rb_ayu

3

Rd_ay

5

Rd_ayg

6

Rd_ayu

2

re_dsr

2

re_pelagic

23

re_pH

22

re_slope

4

Ro_ay

19

Ro_ayg

21

Ro_ayu

32

Rp_ay

3

Rp_ayg

3

Rp_ayu

3

Rs_ay

5

Rs_ayg

11

Rs_ayu

16

Ry_ay

6

Ry_ayg

5

Ry_ayu

13

sd_a_wt

13

sd_wt

14

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI